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Abstract 

We seek evidence of the Yarkovsky effect among Near Earth Asteroids (NEAs) 
by measuring the Yarkovsky- related orbital drift from the orbital fit. To 
prevent the occurrence of unreliable detections we employ a high precision 
dynamical model, including the Newtonian attraction of 16 massive asteroids 
and the planetary relativistic terms, and a suitable astrometric data treat- 
ment. We find 21 NEAs whose orbital fits show a measurable orbital drift 
with a signal to noise ratio (SNR) greater than 3. The best determination is 
for asteroid (101955) 1999 RQ36, resulting in the recovery of one radar ap- 
parition and an orbit improvement by two orders of magnitude. In addition, 
we find 16 cases with a lower SNR that, despite being less reliable, are good 
candidates for becoming stronger detections in the future. In some cases 
it is possible to constrain physical quantities otherwise unknown by means 
of the detected orbital drift. Furthermore, the distribution of the detected 
orbital drifts shows an excess of retrograde rotators that can be connected 
to the delivery mechanism from the most important NEA feeding resonances 
and allows us to infer the distribution for NEAs obliquity. We discuss the 
implications of the Yarkovsky effect for impact predictions. In particular, for 
asteroid (29075) 1950 DA our results favor a retrograde rotation that would 
rule out an impact in 2880. 
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1. Introduction 



It is well known that nongravitational forces should be considered as 
important as collisions and gravi tational perturbatio ns for the overall un- 
derstanding of asteroid evolution ( jBottke et al.l . 120061 ) . The most important 
nongravitational perturbation is the Yarkovsky effect, which is due to radia- 
tive recoil of anisotropic thermal emission and causes asteroids to undergo a 
secular semimajor axis drift da J dt. Typical values of da I dt for sub-kilometer 
NEAs are lO^-lO" 3 au/Myr flVokrouhlickv et al.l . lioooh . 

The Yarkovsky acceleration depends on several ph ysical quantities suc h 
as spin state, size, mass, shape, and ther mal properties (jVokrouhlickyl . Il999l ) . 



Furthermore, iRozitis and Greenl (120121 ) show that surface roughness also 



plays an important role by enhancing the Yarkovsky related semimajor axis 
drift by as much as tens of per cent. Though no complete physical character- 
ization is typically available to compute the Yarkovsky acceleration based on 
a thermophysical model, the orbital drift may be detectable from an astro- 
metric dataset. As a matter of fact, a purely gravitational dynamics could 
result in an unsatisfactory fit to the observational data. This is especially 
true when extremely accurate observations are available, e.g., radar obser- 
vations, or when the observational dataset spans a long time interval thus 
allowing the orbital drift to accumulate and become detectable. 

Until recently, the Yark ovsky effect has bee n measured directly only in 



three cases, (6489) Golevka (ICheslev et al.Ll2003h. (1 52563) 1992 BF flVokrouhlickv et al. 



2008|), and recently for (101955) 1999 RQ 36 flCheslev et all 120121 ). For both 



Golevka and 1999 RQ36 the Yarkovsky perturbation must be included to fit 
accurate radar observations spanning three apparitions. For 1992 BF the 
Yarkovsky effect is needed to link 4 precovery observations of 1953. Fur- 
thermore, in the case of 1999 RQ36 the available physical characterization, 
along with the estimate of the Yarkovsky effect, allows the estimate of the 
ast eroid's bulk d e nsity- 



Nugent et al.l (j2012bl ) find 54 detections of semimajor axis drift by per- 



forming a search for semimajor axis drift among NEAs similar to the one 
presented in this paper. However, there are differences in the observational 
data treatment, in the modeling, and in the selection filters. A description 
of the differences an d a comparison of the results is contained in Sec. 13.21 



Nugent et al.l ( ]2012al ) use WISE-derived geometric albedos and diameters to 



predict orbital drifts for 540 NEAs. Even if none of these objects has an ob- 
servational record that allows one to measure the predicted orbital drift, the 
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authors list upcoming observing opportunities that may reveal the Yarkovsky 
signal. 

The Yarkovsky effect plays an important role for orbital predictions such 
as those concerning Earth impacts. In particular, when an asteroid has an 
exceptionally well constrained orbit, the Ya r kovsk y effect may become the 
principal source of uncertainty. iMilani et al.l (120091 ) show how the size of the 
semimajor axis drift along with its uncertainty modifies impact predictions 
for the next century for 1999 ROj36- The cumulative impact probability is 
approximately 10~ 3 , while a Yarko vsky- free propagation would rule out any 
impact event. IChesley et al.l (120121 ) improve the da/dt estimate by means of 
September 2011 Arecibo radar measurements and find a cumulative impact 
probability approximately 4 x 10~ 4 . Another remarkable c ase is (99942) 



Apoph is. T hough only a margin al da/dt estimate is available, iGiorgini et al 



(l2008h and ICheslev et all (l2009h orove that the occurrence of an impact in 
2036 is decisi yely driven by t he m agnitude of the Yarkovsky effect. On the 
longer term, IGiorgini et al.l ( 120021 ) show that an impact between asteroid 
(29075) 1950 DA and the Earth in 2880 depends on the accelerations arising 
from thermal re-radiation of solar energy absorbed by the asteroid. 



2. Methodology 

2.1. Yarkovsky modeling and determination 

The Yarkovsky effect depends on typically unknown physical quantities. 
As the primary manifestation is a semimajor axis drift, we seek a formulation 
depending on a single parameter to be determined as a result of the orbital 
fit. To bypass the ne e d of physical characterization we used a comet-like 



model (IMarsden et all Il973l ) for transverse acceleration a t = A 2 g(r), where 



g is a suitable function of the heliocentric distance r and A 2 is an unknown 
parameter. To determine A 2 we used a 7-dimensional differential corrector: 
starting from the observational dataset we simultaneously determine a best- 
fitting solution for both the orbital elements and A 2 . 

Once A 2 is determined from the orbital fit we estimate semimajor axis 
drift by means of Gauss' perturbative equations: 



2aVT^7< 



nr 



A 2 g(r) 



(1) 



where a is the semimajor axis, e is the eccentricity and n is the mean motion. 
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By averaging we obtain 



7T 



Jo r una J 



where T is the orbital period and / is the true anomaly. Let us now assume 
g[r) = (r /r) d , where r is a normalizing parameter, e.g., we use r = 1 au. 
In this case the semimajor axis drift is 

S = ^i^^yf(l + ecos/)-V,. (3) 



7m \p 
By Taylor expansion, we have 

/>2tt °° / j _ i\ />2w 

/ (l + ecosf) d - 1 df = J2( . W cos fc M- 
./o fe=0 V ft / ./o 

The odd powers of the cosine average out, so we obtain 

t 2A 2 (l-e 2 ) fr ^ d 



(4) 



n 

where 



(^) (5) 



fc=0 

The ratio 



(7) 

a fc V 2k + 2j \ 2k + 2 J w 

is smaller than 1 for d > and large enough. Therefore, are bounded and 
J(e,d) is convergent for any e < 1. Equation (jZJ) can be used to recursively 
compute ctfc starting from ato — 1. For integer <i the series J is a finite sum 
that can be computed analytically, e.g., J(e, 2) = 1 an d J(e, 3) = 1 + 0.5e^ . 

The proper value of d is not easily determined. From I Vokrouhlickyl ( 119981 ). 
we have 

0.59 ^ o (8) 

9 ^ JJK ' 1 ' JK ' l + e + 0.56 2 V ; 

for the Yarkovsky diurnal component (which is typically dominant) , where 
A is the Bond albedo, is the thermal parameter, 7 is the obliquity, and 



4 



$(r) is the standard radiation force factor, which is inversely proportional 
to the bulk density p, the diameter D, and r 2 . The thermal parameter is 
related to the thermal inertia F by means of the following equation 




where e is the emissivity, a is the Boltzmann's constant, T* is the subsolar 
temperature, and P is the rotation period. In this paper we use d — 2 to 
match the level of absorbed solar radiation. Then, from Eq. (JH]) we have that 

A 2 ~ 4 ^ ~ ^ $(lau)/(9) cos 7 . (10) 
9 

However, as T* oc r~ a5 we have that G oc r L5 , therefore the best value of ti 
depends on the object's thermal properties: 

• for O ^> 1 we obtain / oc r -1 - 5 , which gives d = 3.5; 

• for O 1 we obtain / oc r 15 , which gives d = 0.5. 

These are limit cases, the true d is always going to be between them. As a 
matter of fact, it turns out that most NEAs, whose rotation period is not 
excessively large and whose surface thermal inertia is not excessively small 
or large, have typically values of G near unity or only slight ly larger, and 



we ca n thus expect d value in the range 2-3. As an example, IChesley et al. 



(120121 ) show that for 1999 RQ36 the best match to the Yarkovsky perturbation 
computed by using a linear heat diffusion model is d = 2.75. 

What matters to us is that da/dt does not critically depend on the chosen 
value of d. As an example for asteroid 1999 RQ36 we have that da/dt = 
(-18.99 ±0.10) x 10~ 4 au/Myr for d = 2 and da/dt = (-19.02 ±0.10) x 10~ 4 
au/Myr for d — 3. Another example is Golevka, for which we obtain da/dt = 
(-6.62 ± 0.64) x 10~ 4 au/Myr for d = 2 and da/dt = (-6.87 ± 0.66) x 10" 4 
au/Myr for d = 3. In both cases the difference in da/dt due to the different 
values assumed for d is well within one standard deviation. 

2.2. Dynamical model 

To consistently detect the Yarkovsky effect we need to account for the 
other accelerations down to the same order of magnitude. For a sub-kilometer 
NEA, typical values of a t range from 10~ 15 to 10~ 13 au/d 2 . 
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Table 1: Gravitational parameters of perturb i ng ast eroids. The masses of Ceres, Pal- 
las, and Vesta are from IStandish and Hellingd (l989), the ones of Euphrosyne and Her- 
culina are from the Institute of Applied Astronomy of RAS, St. Petersburg, Russia 
(http ://www.ipa.nw. r u/PA GE/DEPFUND/LSBSS/cngmasses. htm) the mass of Juno is 
from iKonopliv et al.l (|201ll ). all the others are from lBaer et al.l (|201ll ). 



Astproid 


Crin fkm^/s 2 l 

V — X 1.1.1. IVlll / k.j 


ill I ,PfPQ 


uo.^uu 


(9^ Pallas 


L l ±.0\J\J 


(4) Vesta 


17.800 


(10) Hygea 


6.0250 


(29) Amphitrite 


1.3271 


(511) Davida 


3.9548 


(65) Cybele 


1.0086 


(9) Metis 


1.3669 


(15) Eunomia 


2.2295 


(31) Euphrosyne 


1.1280 


(52) Europa 


1.2952 


(704) Interamnia 


4.7510 


(16) Psyche 


1.7120 


(3) Juno 


1.9774 


(532) Herculina 


1.5262 


(87) Sylvia 


1.3138 



Our N-body model includes the Newtonian accelerations of the Sun, eight 
planets, the Moon, and Plu to that are based on JPL's DE405 planetary 
ephemerides (IStandishl . |2000| ). Furthermore, we added the contribution of 16 
massive asteroids, as listed in Table [TJ 

We used a full relativistic force model including the contribution of the 
Sun, the planets, and the Moon. Namely, we used the Ei nstein - Infeld - 
Hoffman (EIH) app roximation as describ ed in iMoyerl (120031 ) or IWilll (119931 ). 
As already noted in IChesley et al.l (120121 ). the relativistic term of the Earth 
should not be neglected because of significant short range effects during Earth 
approaches that NEAs may experience. For asteroids with a large perihelion 
distance such as Golevka also Jupiter's term could be relevant. Figure [1] 
compares the relativistic accelerations of the Earth and Jupiter as they for- 
mally appear in the EIH equations of motion to the transverse Yarkovsky 
acceleration acting on Golevka. 
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Figure 1: Relativistic accelerations of the Earth and Jupiter as they formally appear in 
the EIH equations of motion compared to the transverse Yarkovsky acceleration acting on 
Golevka. 



2.3. Observational error model 

The successful detection of the Yarkovsky effect as a result of the orbital 
fit strongly depends on the quality of the observations involved. In particu- 
lar, the availability of radar data is often decisive due to the superior relative 
accuracy of radar data with respect to optical ones. Moreover, radar mea- 
surements are orthogonal to optical observations: range and range rate vs. 
angular position in the sky. 

Since the Yarkovsky effect acts as a secular drift on semimajor axis we 
have a quadratic effect in mean anomaly: the longer the time span the 
stronger the signal. However, the presence of biases in historical data and 
unrealistic weighting of observations may lead to inaccurate results. To deal 
wi th this problem we a pplied the debiasing and weighting scheme described 



in 



Chesley et al.l (j2010h . This scheme is a valid error model for CCD ob- 
servations, while for pre-CCD data the lack of star catalog information and 
the very uneven quality of the observations represents a critical problem. In 
these cases the occurrence of unrealistic nominal values for Yarkovsky model 
parameters presumably point to bad astrometric treatments and have to be 
rejected. 

To prevent outliers from sp oiling orbital fi t s, we applied the outlier rejec- 



tion procedure as described in ICarpino et al.l (120031 ). 



Besides the astrometric treatment described above, in the following cases 
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we applied an ad hoc observation weighting: 



1999 RQ36: as already explained by lChesley et al.l ( 120121 ). in some cases 



there are batches containing an excess of observations from a single 
observatory in a single night. To reduce the effect of these batches to a 
preferred contribution of 5 observations per night, we relaxed the weight 
by a factor a/ iV/5, where N is the number of observations contained 
in the batch. 

1992 BF: as the precovery observation s of 1953 have been carefully 



remeasured in lVokrouhlicky et al.l (120081 ) . these observations were given 



a weight 0.5" in right ascension and 1" in declination. 

• Apollo: the large dataset available for Apollo contains observations go- 
ing back to 1930. Many pre-CCD era observation batches show unusu- 
ally high residuals, especially during close Earth approches. To lower 
the effect of non-CCD observation, we used weights 10" for observations 
from 1930 to 1950 and 3" from 1950 to 1990. 

• 1989 ML: the discovery apparition contains observations from Palomar 
Mountain showing large residuals whether or not the Yarkovsky per- 
turbation is included in the model. Even if this apparition increases 
the observed arc by three years only, we felt it safer to weight the 
corresponding observations 3" . 

2.4- Treatment of precovery observations 

There are a few cases where the Yarkovsky signal is mainly contained in 
few isolated precovery observations. This is the case of the already mentioned 
asteroid 1992 BF, which has 4 isolated observations in 1953 from Palomar 
Mountain DSS. Other cases are 

• 1999 FK 2 i, which has 6 isolated observations in 1971 from Palomar 
Mountain; 

• 2001 MQ 3 , which has 4 isolated observations in 1951 from Palomar 
Mountain DSS; 

• 1989 UQ, which has 3 isolated observations in 1954 from Palomar 
Mountain; 
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1999 FK 21 



2001 MQ 3 




Figure 2: Observed position, average observed position, and predicted postfit position with 
the Yarkovsky perturbation in the dynamics with respect to the predictions without the 
Yarkovsky perturbation in the dynamics. The uncertainty ellipses for the mean observed 
position and the predicted positions correspond to the 3-tr level. The dashed lines represent 
the predicted motion in the sky of the asteroid in the precovery apparition. Circles are the 
predicted positions, crosses are the measured positions, stars are the mean of the measured 
positions, while diamonds are the predicted positions when the Yarkovsky perturbation is 
included in the dynamics. The origin is arbitrarily set to the non- Yarkovsky prediction. 



1991 VE, which has 4 isolated observations in 1954 from Palomar 
Mountain DSS. 

For all these cases it would be desirable to remeasure th e precovery obser- 

" (120081 ) 



vations as done for 1992 BF in iVokrouhlicky et al 



where precovery 

observations were corrected by an amount up to 3.1" . For this reason we con- 
servatively gave weights 3" to the precovery observations of the four asteroids 
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above. 

Besides the conservative weighting, we ruled out clock error as a possible 
cause of the Yarkovsky signal. Figure |2] shows the scenario for the four men- 
tioned asteroids during the precovery apparition. We can see that it is not 
possible to match the observations by translating the non- Yarkovsky uncer- 
tainty ellipse on the along track direction. The Yarkovsky solution produces 
a shift in the weak direction that give a better match to the observations, in 
particular when we take the average of the observed positions. 

2.5. Filtering spurious results 

To assess the reliability of the Yarkovsky detections we computed an 
expected value for A 2 starting from the 1999 RQ36 case, which is the strongest 
and most reliable detection, by scaling according to ( ITU]) 

{M) exp = (A 2 ) RQm ^g* . (ii) 

For diameter D we used either the known value when available or an assumed 
value compu ted from the absolute ma gnitude H according to the following 



relationship (jPravec and Harris! . 120071 ): 



IQ-0.2H 

D = 1329 km x (12) 

VPv 



where py is the albedo, assume d to be 0.154 if unkn own, in agreement with 



the Palermo scale computation (IChesley et al.l . 12002] ). 



Some physical properties of 1999 RQ36 maximize A 2 (IChesley et al.l . 12012] ). 
In particular 7 is nearly 180°, A is 0.01, and p is low (0.96 g/cm 3 ). On 
the other hand G = 4.33 for which /(0) — 0.15 while the maximum is 
~ 0.21. For these reasons we selected those Yarkovsky detections for which 
S = \A2/(A2)exp\ was smaller than 1.5. The selected threshold allows some 
tolerance as we are scaling only by D without accounting for other quanti- 
ties such as bulk density, thermal properties, obliquity, spin rate, and surface 
roughness. 

A high SNR threshold is likely to produce robust detections with respect 
to the astrometric data treatment. With lower SNR the sensitivity to the 
observation error model increases and detections become less reliable. We 
decided that 3 was a sensible choice for minimum SNR, even if we analyze 
detections with smaller SNR in Sec 
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3. Results 

We applied our 7- dimensional differential corrector to determine the pa- 
rameter A 2 and the corresponding da/dt for all known NEAs. After applying 
the filters discussed in Sec. 12.51 we obtain 21 Yarkovsky detections that we 
consider reliable (Table [2]). The reported uncertainties are marginal, i.e., 
they fully take into account the correlation between A 2 (and thus da/dt) and 
the orbital elements. 

We cross-checked these detections by using two independent software 
suites: the JPL Comet and Asteroid Orbit Determination Package and Orb- 
Fit (http://adams.dm.unipi.it/orbfit/ 0: in all cases we found agreement at 
better than the 1-a level. 



3.1. 2009 BD 

Asteroid 2009 BD is very small and to fit its observational dataset it is 
necessary to include solar radiation pressure ( iMicheli et all 120121 ) . Therefore, 
we also included in the model a radial acceleration a r = Ai/r 2 . Along with 
the tabulated value of A 2 , we obtained A\ = (62.05 ± 8.85) x 10~ 12 au/d 2 . 
This results in an area to mass ratio A/M = (2.72 ± 0.39) x 10~ 4 m 3 /kg, 



which is consistent at the 1-a level with the value reported by iMicheli et al 
f l2012h . i.e., (2.97 ±0.33) x 10~ 4 m 3 /kg. 



3.2. Comparison with other published results 

The first three objects of Tableware the already known cases of Golevka, 
1992 BF, and 1999 RQ 36 . Whi le for 1999 RQ^ th ere is a perfect match 



between our result and the one i nlCheslev et al.l ( 20121) . f or Golevka and 1992 
BF th e values are different from lChesley et al.l ( 120031 ) and lVokrouhlicky et al. 
( 120081 ) . respectively. However, this can be easily explained by the availability 
of new astronomy and the fact that the present paper adopted the debiasing 
and weighting scheme by IChesley et al.l ( 2010f ). which was not available at 
the time of the earlier pub lications. 



As already mentioned, iNugent et al.l ( I2012bl ) performed a search similar 
to the one presented in this paper and found 54 NEAs with a measurable 
semimajor axis drift. The main differences are the following: 



OrbFit was used in the development version 4.3, currently in beta-testing. 
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• They selected only numbered objects, while we included all known 
NEAs. 

• Their observation dataset was slightly different as they used observa- 
tions until 2012 January 31, while we have data until 2012 October 31. 
This difference does not really matter for optical data, but it does for 
radar data, e.g., for Aten and Toro. Moreover, they did not use single 
apparition radar, while we did as we think they represent an important 
constraint. 

• They solved for constant da/dt while we used constant A2 and then 
convert to da/dt. These techniques are equivalent when the semimajor 
axis and eccentricity are constant, but there could be differences as we 
cannot assume da/dt constant for objects experiencing deep planetary 
close approaches. 

• They searched for the best-fit da/dt by means of the golden section 
algorithm, i.e., they computed the RMS of the residuals correspond- 
ing to the best fitting orbital elements for fixed values of da/dt, then 
the minimum is obtained by interpolation. In this paper we used a 
full 7-dimensional differential corrector. The two methods should be 
equivalent. 

• They used 1 as lower bound for SNR, while we use 3 that gives detec- 
tions more robust against changes in the observation weighting. Also, 
they used the "sensitivity" parameter, i.e., a metric to measure the 
sensitivity of a dataset to the presence of a semimajor axis drift. We 
did not use such a metric as we think that an SNR> 3 is already a 
good metric in that respect. 

• We kept only those objects for which the measured orbital drift can be 
related to the Yarkovsky perturbation presuming that inconsistencies 
stem from astrometric errors, while they also considered the possibility 
of other nongravitational effects such as a loss of mass. 

Table [2] contains a compariso n between our orbital drifts and the ones 



reported bv iNugent et all (|2012bJ). 2009 BD, 1999 MN, 1999 FA, and 2005 
ES70 are present only in our list as they are not numbered, while 2001 MQ3, 
1989 UQ, and 1991 VE are eliminated by their filters. It is worth pointing 
out that also 1999 RQ36, 1992 BF, Golevka, and YORP have been filtered out 
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by lNugent et al.l (j2012bl ) criteria, even though t h ey rep ort the corresponding 
detections for a compariso n with IChesley et al.l (120081) . 

Among the cases that iNugent et al.l (j2012bl ) report with SNR > 3 we 
neglected the following three: 

• (1036) Ganymed for which we found A 2 = (-16.54±4.35) x 10~ 15 au/d 2 , 
corresponding to da/dt = (—6.06 ±1.59) x 10 -4 au/Myr. However, the 
nominal A 2 is 28 times larger than (A 2 ) exp , so we marked this detection 
as spurious. As Ganymed observations go back to 1924, this unreliable 
detection might be due to bad quality astrometry. 

• (4197) 1982 TA for which we used the radar apparition of 1996, which 
reduced the SNR below 1. For this object we found A 2 = (5.61 ± 
14.26) x 10~ 15 au/d 2 , co rresponding to da/dt = (3.88 ± 9.86) x 10~ 4 
au/Myr. For comparison INugent et al.l (j2012bl ) report da/dt = (30.9 ± 
9.2) x 10~ 4 au/Myr. 

• (154330) 2002 VX 4 for which we found A 2 = 102.36 ± 36.34) x 10~ 15 
au/d 2 , corresponding to da/dt = (43.10±15.25) x 10~ 4 au/Myr. Again, 
the value of A 2 was ~ 4 times larger than (A 2 ) exp , so also this detection 
was marked as spurious. 

Besides the differences outlined above, there is an overall agreement for 
the common cases. As a matter of fact by computing the a of the difference, 
i.e., cr 2 = o\ + u 2 , there are only two cases that are not consistent at the 1-a 



level: 



1992 BF, for which we used remeasured observations (jVokrouhlicky et al. . 
20081 ) and changed the weights (see Sec. 12. 3p for the 1953 apparition. 
The remeasurements are not available in the MPC database. 

Apollo, for which we applied a suitable manual weighting as described 
in Sec. 12.31 



3.3. Lower SNR and other less reliable detections 

Table [3] contains detections that we rate as nonspurious on the basis of 
the S ratio between expected and measured value, down to SNR = 2 plus 
the following remarkable cases: 
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Figure 3: Same of Fig. [2] for asteroid 1999 FA. The uncertainty ellipses for predicted posi- 
tions correspond to 3-cr levels, while the uncertainty for the observaed position correspond 
to 3". 
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Figure 4: Relative displacement of A2 with respect to the nominal value as a function of 
the number of perturbers for 1999 RQ36 and Toutatis. 

• 1999 FA, which has 1 isolated observation in 1978 from Siding Spring 
Observatory. By performing the same analysis of Sec. 12.41 (see Fig. [3]) 
we see no clear improvement due to Yarkovsky as the observation is 
3.5-4 a away from the prediction either way. In this case a clock error 
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may bring the Yarkovsky solution close to the observation. As this 
detection depends on a lone, isolated observation we would rather be 
cautious and consider this detection reliable only if the 1978 observation 
is remeasured. 

Eros, which looks like a relia ble detection as S = 0.75. However, the 
obliquity is known to be 89° ( Yeomans et al. . 2000l ) and then we enter 



the regime where the seasonal component of the Yarkovsky effect is 
dominant. As the seasonal componen t is typ ically 10 times smaller then 



the diurnal one (jVokrouhlicky et al.l . 120001 ) we mark this detection as 



spurious. This points to possible bad astrometric treatment, especially 
for historical observations dating back to 1893. 

Toutatis, which enters the Main Belt region because of the 4.12 au 
aphelion and the low inclination (0.44°). Therefore, it is important 
to account for the uncertainty in the masses of the perturbing aster- 
oids. By taking into account the uncertainty of the perturbing asteroid 
masses the actual uncertainty in A 2 increases by 11% with a commen- 
surable drop in SNR. Figure H] shows the evolution of A 2 as a function of 
the number of perturbers for asteroids 1999 RQ36 and Toutatis. While 
for 1999 RQ36 (aphelion 1.36 au) we reach convergence with just four 
perturbers, for Toutatis we have a quite irregular behavior suggesting 
that we may need to include more perturbers. 

1994 XLi, whose observations in 1994 from Siding Spring Observatory 
show high residual so we relaxed weights to 3" . 

2005 QC 5 and 2000 NL 10 have been included despite the low SNR. Sim- 
ilarly to the cases described in Sec. \2A\ we applied weights 3" to pre- 
covery observations and this data treatment weakened the Yarkovsky 
detection. However, we think that remeasuring the precovery observa- 
tions may lead to reliable Yarkovsky detections for these objects. 

1950 DA, for which the Yarkovsky effect plays an important role for 
impact predictions, e.g., see Sec. 14.41 

Icarus, for which the low values of both the nominal value and the un- 
certainty put a strong contraint on A 2 . In particular , it is worth point- 



i ng ou t that the constraint on da/dt is consistent with lVokrouhlicky et al 
f l2000l . Fig. 5) where \da/dt\ < 3 x 10~ 4 au/Myr. 
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There are other objects with an even lower SNR for which the Yarkovsky 
signal might be revealed if precovery observations were remeasured: (11284) 
Belenus, (66400) 1999 LT 7 , (4688) 1980 WF, (67399) 2000 PJ 6 , (267759) 
2003 MC 7 , and (88710) 2001 SL 9 . 

Though these detections have to be considered less reliable, some of them 
may be good candidates for becoming stronger detections in the fu t ure if high 
quality astrometry is obtained, e.g., by radar or GAIA ( iMignardl . 120031 ). 



4. Discussion 



4-1. Connection with NEA feeding mechanisms 

The diurnal Yarkovsky effect produc es a semimajor axis drift propor- 



tional to cos 7 ( jVokrouhlicky et all 120001 ) . As the diurnal term is typically 



the dominant one, the sign of da/dt can be related to the asteroid spin ori- 
entation, i.e., a negative da/dt corresponds to a retrograde rotator while a 
positive da/dt corresponds to a prograde rotator. This conclusion is sup- 
ported by the eight known obliquities for the asteroids in the sample that 
are listed in Table [2j in all cases the spin axis obliquity is consistent with 
the sign of da/dt. 

We can now use this interpretation and our solution for the Yarkovsky 
semimajor axis drift values for NEAs in the following way. Table H] contains 
four prograde rotators and seventeen retrograde rotators. This excess of 
retrograde rotators can b e explained by the n ature of resonance feeding into 



the inner Solar System ( IBottke et al.l . 120021 ) . Most of the primary NEA 
source regions (e.g., the 3:1 resonance, JFCs, Outer Belt, etc.) allow main 
belt asteroids to enter by drifting either inwards or outwards, but the v§ 
resonance is at the inner edge of the main belt and so ast eroids can generally 
enter only by inwards drift, i.e., with retrograde r otation. IBottke et al.l (120021 ) 
report that 37% of NEAs arrive via v§ resonance. La Spina et al.l ( 2004 ) point 
out that this implies 37% of NEAs have retrograde spin (via z/ 6 ), plus half of 
the complement (via other pat hways). Thus, the retr ograde fraction should 
be 0.37+0.5x0.63 = 0.69, while lLa Spina et al.l Jiool report 67% retrograde 
for their sample, which is dominated by large NEAs. 

Table [2] contains 81% retrograde rotators, which is larger than 69% and 
then appears to be inconsistent with the theory. However, the asteroid sam- 
ple of Table |5] is based on measured Yarkovsky mobility, and so is dominated 
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Figure 5: Left panel: distribution of A 2 and its uncertainty. Right side of plot is prograde 
rotation and left side is retrograde rotation. Right panel: running box mean in SNR for 
the fraction of asteroids with indication of retrograde spin. Cases with SNR > 2 show a 
consistent 4:1 retrograde ratio, cases with little or no signal are split 50-50. 



by small PHAs. For this populatio n, 55% are expected to arrive from v§ 



according to the lBottke et al.l (120021 ) NEO model (Bottke, personal commu- 
nication), a much greater fraction than for NEAs as a whole. The same 
arithmetic used for NEAs gives 0.55+0.5x0.45=0.78, which matches very 
well the 81% ± 8% value of Table H 

To assess the behavior of the fraction of retrograde rotators as a function 
of the SNR, we took all of the objects with S < 1.5. The left panel of Fig. [5] 
shows the distribution of A 2 and its uncertainty. The excess of retrograde 
rotators is clearly visible for SNR > 3 and also between 2 and 3, where 
small PHAs dominate. For lower SNR we have a more uniform distribution. 
The right panel of Fig. [5] is a running mean of the fraction of retrograde 
rotators as a function of the SNR. For SNR < 1 we are in a noise dominated 
regime for which we have a rough 50% fraction of retrograde rotators, for 
1 < SNR < 2 we have a transition from noise-dominated to signal dominated, 
and for SNR > 2 we have a signal dominated regime with around 80% 
retrograde rotators. 

We can also try to use the detected values to infer the obliquity distribu- 
tion. From Eq. (TTUT) we have that 

A 2 oc (13) 

and so we can estimate 7 by using the either known or estimated by Eq. ffT21 
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Figure 6: Inferred obliquity distribution for objects with SNR > 1. 



diameter and using 1999 RQ36- By taking all the detections with SNR > 1 
and assuming a fixed p = 1500 kg/m 3 we obtain the distribution of Fig. El 
where the cases with | 00S7I > 1 have been placed in the extreme obliquity 
bins. Despite the low number of bins, we can see the excess of retrograde 
rotators and the abund ance of objects with an ex t reme obliquity, as expected 
from the YORP effect ( jCapek and Vokrouhlickyi . I2004T ) . While this distribu- 
tion should be considered only approximate due to the numerous assumptions 
(e.g., neglecting dependence on bulk density, shape and thermal properties) 
we consider it to be a significant improvement over what is otherwise known. 
However, we find it interesting that it app ears to be cons i stent with the 
observed obliquity distribution of the NEAs (ILa Spina et al.l . 12004 ). 



4-2. Spurious detections 

Our search for Yarkovsky signal produced a large number of spurious de- 
tections, i.e., semimajor axis drifts far larger than the Yarkovsky effect would 
cause. Figure [7] contains the histograms of S for different SNR intervals. For 
SNR > 3 we have 67% spurious detections, for 2 < SNR < 3 we have 88%. 

It it worth trying to understand the reason of these spurious solutions. 
We think there are two possible causes: 



Bad astrometry treatment: as discussed in Sec. 12.31 and Sec. \2A\ non- 
CCD observations may contain errors that are difficult to model. If an 
ad hoc weighting is not used we may have misleading results. Indeed, 
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Figure 7: Histogram of S for different intervals of SNR. The vertical lines correspond to 
the selected threshold 1.5. 



spot checking of such cases generally confirmed isolated astrometry as 
source of spurious detections. 

• Incompleteness or inconsistency in the dynamical model: the formula- 
tion proposed in Sec. 12. H is a simplified model of the Yarkovsky pertur- 
bation that might be poor in some cases. A more sophisticated formu- 
lation would require a rather complete physical characterization that 
is typically unavailable and thus cannot be used for a comprehensive 
search as done in this paper. Moreover, as discussed in Sec. I3.3l for Tou- 
tatis, we may need to include more perturbing asteroids (and the un- 
certainty in their masses) in the model. Finally, we cannot rule out the 
possibility of nongravitational perturbations different from Yarkovsky. 



4-3. Constraining physical quantities 

The results reported in Table [2] can be used to constrain physical quan- 
tities. When A, D, and 7 are known Eq. ffTUl) provides a simple relationship 
between p and 0. This relationship can be easily translated into a rela- 
tionship between p and the thermal inertia V by means of Eq. (Q. As a 
benchmark of these technique we c an use asteroid 1999 RQ36 (Fig. [8]), for 
which the known values of T and p ( jChesley et al.l . |2012| ) match the plotted 
contraint. Figure M shows the possible values of p as a function of T for as- 
teroids Golevka, Apollo, Ra-Shalom, Toro, YORP, and Geographos. For the 
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Figure 8: Density as a function of thermal inertia for asteroid 1999 RQ36- The three lines 
corresponding to the nominal value and the 1-a levels are extremely narrow because of the 
high SNR. The circ le corresponding to t he nominal observed value of T and thus inferred 



nominal value of p (jCheslev et all 120 



H) matches the plo tted curves. The enclosed region 



corresponds to the 1-a interval for T ([Emery et al 
the 1-a limits of p. 



2012) and indicates how it maps onto 



latter two objects we assumed slope parameter G = 0.15. Our findings are 
consistent with the taxonomic type. For instance, Golevka, Apollo, Toro, 
YORP, and Geographos are S/Q-type asteroids with an expected density 
between 2000 and 3000 kg/m 3 , while Ra-Shalom is a C-type so we expect a 
bulk density from 500 and 1500 kg/m 3 . Figure [9] suggests that Golevka has 
thermal inertia 150 < T < 500 J m~ 2 s~ ' 5 K - 1 and Apollo has a rather large 



thermal inertia 400 < T < 1000. According to lDelbo et al.l (120031 ) Ra-Shalom 
has an unusually high thermal inertia. In fact, by taking the right side of 
the plotted rectangle we obtain a density closer to 1000 kg/m 3 , similar to 
the one 1999 RQ36, which belongs to a similar taxonomic class. For Toro, 



Muellerl (120121 ) reports a thermal inertia 200 < T < 1200 J m" 2 s~ 5 K" 



but likely lower, which would result in a bulk density between 2000 and 4000 
kg/m 3 . 

4-4- Implications for impact predictions 

The question is how much the predictions of possible asteroid impacts on 
Earth can be affected by the Yarkovsky effect. There are already three well 
known examples, namely the asteroids Apophis, 1999 RQ36, and 1950 DA 
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for which Yarkovsky perturbations are relevant and need to be accounted for 
in the impact risk assessment. 

For Apophis and 1999 RQ36 this is due to the presence of a strongly 
scattering planetary close approach between now and the epochs of the pos- 
sible impacts. These encounters transform a very well determi ned orbit into 



a po orly known one for which chaotic effects are dominant (IMilani et al 



20091 ). Apophis will have a close approach in April 2029 with minimum 
distance of ~ 38000 km from the geocenter. As a consequence the orbital 
uncertainty will increase by a factor > 40000. 1999 RQ36 will have a close 
approach to Earth in 2135 with nominal minimum distance about the same 
as the distance to the Moon, with an increase in uncertainty by a factor 
~500. The minimum possible distance for this close approach is three times 
smaller and would re sult in an increase of the uncertainty by a factor ~1500 
(jChesley et all |2012| ). In both cases the Yarkovsky effect is large enough 
to shift the position at the scattering close a pproach by an amount much 
larger than the distance between the keyholes ( Chodasl . 1999 ) corresponding 
to impacts in later years (2036, 2037, 2068 for Apophis; 2175, 2180, 2196 
for 1999 RQ36). Thus, the occurrence of these later impacts is determined 
by the Yarkovsky perturbation in the years between now and the scattering 
encounter. For 1950 DA the influence of the Yarkovsky effect for the possible 
impact in 2880 is due to long time int erval preceding the i mpact that allows 
the orbital displacement accumulate (jGiorgini et al.L 120021 ) . 

Currently, 1999 RQ36 is the case with the best determined Yarkovsky ef- 
fect (SNR ~ 200), while Apophis and 1950 DA have only a marginal detection 
(SNR < 1 and SNR ~ 1.4, respectively). Therefo re, the impact monit oring; 
for 1999 RQ36 fully takes into account Yarkovsky ( Chesley et al.L 2012 ). On 
the contrary, the current estimate of impact probabilities for Apophis is based 
on a Montecarlo model of Yark ovsky based on a pri ori knowledge of the sta- 
tistical properties of this effect (jChesley et al.L 120091 ) . 1950 DA could benefit 
from a similar approach. The knowledge of the Yarkovsky perturbation for 
Apophis is expected to increase dramatically due to observations possible 
during the January 2013 close approach. 

We investigated the possibility that our identification of asteroids with 
measurable Yarkovsky produces new cases such as the two above, that is 
of impact monitoring affected by Yarkovsky. The answer to this question 
is negative, in that the intersection between the current list of NEA with 
possible impacts on Earth (337 according to NEODyS, 404 according to 
Sentry) and the list with detected Yarkovsky contains only 1999 RCj36- 
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However, this conclusion depends on the fact that our monitoring of 
possible future impacts is done for only about one century (currently 90 
years for NEODyS, 100 for Sentry). 1999 RQ36 was a special case, re- 
lated to an intensified effort for the OSIRIS-REx mission target 1999 RQ36 



(http:/ /osiris-rex.lpl.arizona.edu). If this time span were generally increased 
to 150-200 years, there co uld well be o t her c ases similar to 1999 RQ36. 



For asteroid 1950 DA iBusch et al.l (120071 ) report a density around 3000 
km/m 3 and two possible solutions for pole orientation and effective diameter: 

1. 7 = 24.47° and D = 1.16 km; 

2. 7 = 167.72° and D = 1.30 km. 

By scaling from the 1999 RQ36 case, we obtain A 2 = 7.01 au/d 2 for the 
direct solution and A 2 = —5.83 au/d 2 for the retrograde solution. Even if we 
found a low SNR detection, our result strongly favors the retrograde solution, 
which is at 0.6 a, than the direct solution, which is more than 3 a away. As 



a consequence, the 2880 impact would be ruled out (jGiorgini et al.l . 120021 ). 



5. Conclusions 

In this paper we developed a 1-parameter formulation that models the 
transverse component of the Yarkovsky effect. Despite being simple, this 
formulation captures the essence the Yarkovsky effect and also allows us to 
solve for the introduced parameter as a part of the orbital fit. 

The small magnitude of the Yarkovsky acceleration made it necessary 
to employ a high precision dynamical model including the effects of sixteen 
perturbing asteroids and planetary relativistic terms. 

As the solution of the orbital fit depends on the obs ervational data treat - 
ment, we applied the debiasing and weighting scheme by lChesley et al. 



Moreover, we resorted to some manual weighting in cases where an incorrect 
treatment non-CCD observations could have led to unreliable results. 

After filtering out those objects with an unreasonably high orbital drift, 
we obtained 21 robust detections with a signal to noise ratio grater than 
3. By connecting the sign of the orbital drift to the spin orientation, we 
found a 4:1 fraction of retrograde rotators. We discussed how this excess of 
retrograde rotators can be related to the delivery of NEAs to the inner Solar 
System. 

For a few asteroids we were able to use the measured orbital drift along 
with the known physical characterization to constrain unknown physical 
quantities such as bulk density and thermal inertia. 
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We also discussed the implication of the Yarkovsky effect and pointed 
out how this perturbation can be relevant either when there is a scattering 
close approach before the possible impact or when the impact is far enough 
in time. 

Acknowledgments 

DF was supported for this research in part by an appointment to the 
NASA Postdoctoral Program at the Jet Propulsion Laboratory, California 
Institute of Technology, administered by Oak Ridge Associated Universities 
through a contract with NASA, and in part by ESA/ESTEC-SpaceDyS Ser- 
vice Level Agreement SSA-NEO-ESA-SLA2-002_NEODYS. 

SC conducted this research at the Jet Propulsion Laboratory, California 
Institute of Technology, under a contract with the National Aeronautics and 
Space Administration. 

The work of DV was partially supported by the Czech Grant Agency 
(grant 205/08/0064) and Research Program MSM0021620860 of the Czech 
Ministry of Education. 

References 

Baer, J., Chesley, S. R., Matson, R. D., May 2011. Astrometric Masses of 
26 Asteroids and Observations on Asteroid Porosity. The Astronomical 
Journal 141, 143. 

Bottke, W. F., Morbidelli, A., Jedicke, R., Petit, J.-M., Levison, H. F., 
Michel, P., Metcalfe, T. S., Apr. 2002. Debiased Orbital and Absolute 
Magnitude Distribution of the Near-Earth Objects. Icarus 156, 399-433. 

Bottke, Jr., W. F., Vokrouhlicky, D., Rubincam, D. P., Nesvorny, D., May 
2006. The Yarkovsky and Yorp Effects: Implications for Asteroid Dynam- 
ics. Annual Review of Earth and Planetary Sciences 34, 157-191. 

Bowell, E., Hapke, B., Domingue, D., Lumme, K., Peltoniemi, J., Harris, 
A. W., 1989. Application of photometric models to asteroids. In: Binzel, 
R. P., Gehrels, T., Matthews, M. S. (Eds.), Asteroids II. pp. 524-556. 

Brozovic, M., Ostro, S. J., Benner, L. A. M., Giorgini, J. D., Jurgens, R. F., 
Rose, R., Nolan, M. C, Hine, A. A., Magri, C, Scheeres, D. J., Margot, 



23 



J.-L., May 2009. Radar observations and a physical model of Asteroid 4660 
Nereus, a prime space mission target. Icarus 201, 153-166. 

Busch, M. W., Giorgini, J. D., Ostro, S. J., Benner, L. A. M., Jurgens, R. F., 
Rose, R., Hicks, M. D., Pravec, P., Kusnirak, P., Ireland, M. J., Scheeres, 
D. J., Broschart, S. B., Magri, C, Nolan, M. C, Hine, A. A., Margot, 
J.-L., Oct. 2007. Physical modeling of near-Earth Asteroid (29075) 1950 
DA. Icarus 190, 608-621. 

Capek, D., Vokrouhlicky, D., Dec. 2004. The YORP effect with finite thermal 
conductivity. Icarus 172, 526-536. 

Carpino, M., Milani, A., Chesley, S. R., Dec. 2003. Error statistics of asteroid 
optical astrometric observations. Icarus 166, 248-270. 

Chesley, S. R., Baer, J., Monet, D. G., Nov. 2010. Treatment of star catalog 
biases in asteroid astrometric observations. Icarus 210, 158-181. 

Chesley, S. R., Chodas, P. W., Milani, A., Valsecchi, G. B., Yeomans, D. K., 
Oct. 2002. Quantifying the Risk Posed by Potential Earth Impacts. Icarus 
159, 423-432. 

Chesley, S. R., Milani, A., Tholen, D., Bernardi, F., Chodas, P., Micheli, M., 
Sep. 2009. An Updated Assessment Of The Impact Threat From 99942 
Apophis. In: AAS/Division for Planetary Sciences Meeting Abstracts #41. 
Vol. 41 of AAS/Division for Planetary Sciences Meeting Abstracts, p. 
#43.06. 

Chesley, S. R., Nolan, M. C, Farnocchia, D., Milani, A., Emery, J., Vokrouh- 
licky, D., Lauretta, D. S., Taylor, P. A., Benner, L. A. M., Giorgini, J. D., 
Brozovic, M., Busch, M. W., Margot, J.-L., Howell, E. S., Naidu, S. P., 
Valsecchi, G. B., Bernardi, F., May 2012. The Trajectory Dynamics of 
Near-Earth Asteroid 101955 (1999 RQ36). LPI Contributions 1667, 6470. 

Chesley, S. R., Ostro, S. J., Vokrouhlicky, D., Capek, D., Giorgini, J. D., 
Nolan, M. C, Margot, J.-L., Hine, A. A., Benner, L. A. M., Chamber- 
lin, A. B., Dec. 2003. Direct Detection of the Yarkovsky Effect by Radar 
Ranging to Asteroid 6489 Golevka. Science 302, 1739-1742. 

Chesley, S. R., Vokrouhlicky, D., Ostro, S. J., Benner, L. A. M., Margot, J.- 
L., Matson, R. L., Nolan, M. C, Shepard, M. K., 2008. Direct Estimation 



24 



of Yarkovsky Accelerations on Near-Earth Asteroids. LPI Contributions 
1405, 8330. 

Chodas, P. W., 1999. Orbit uncertainties, keyholes, and collision probabili- 
ties. In: Bulletin of the American Astronomical Society. Vol. 31 of Bulletin 
of the American Astronomical Society, p. 1117. 

Delbo, M., DelPOro, A., Harris, A. W., Mottola, S., Mueller, M., Sep. 2007. 
Thermal inertia of near-Earth asteroids and implications for the magnitude 
of the Yarkovsky effect. Icarus 190, 236-249. 

Delbo, M., Harris, A. W., Binzel, R. P., Pravec, P., Davies, J. K., Nov. 2003. 
Keck observations of near-Earth asteroids in the thermal infrared. Icarus 
166, 116-130. 

Drummond, J. D., Wisniewski, W. Z., Feb. 1990. The rotational poles and 
shapes of 1580 Betulia and 3908 (1980PA) from one apparition. Icarus 83, 
349-359. 

Durech, J., Vokrouhlicky, D., Baransky, A. R., Breiter, S., Burkhonov, O. A., 
Cooney, W., Fuller, V., Gaftonyuk, N. M., Gross, J., Inasaridze, R. Y., 
Kaasalainen, M., Krugly, Y. N., Kvaratshelia, O. I., Litvinenko, E. A., Ma- 
comber, B., Marchis, F., Molotov, I. E., Oey, J., Polishook, D., Pollock, 
J., Pravec, P., Sarneczky, K., Shevchenko, V. G., Slyusarev, I., Stephens, 
R., Szabo, G., Terrell, D., Vachier, F., Vanderplate, Z., Viikinkoski, M., 
Warner, B. D., Nov. 2012. Analysis of the rotation period of asteroids 
(1865) Cerberus, (2100) Ra-Shalom, and (3103) Eger - search for the 
YORP effect. Astronomy & Astrophysics 547, A10. 

Durech, J., Vokrouhlicky, D., Kaasalainen, M., Higgins, D., Krugly, Y. N., 
Gaftonyuk, N. M., Shevchenko, V. G., Chiorny, V. G., Hamanowa, H., 
Hamanowa, H., Reddy, V., Dyvig, R. R., Oct. 2008a. Detection of the 
YORP effect in asteroid (1620) Geographos. Astronomy & Astrophysics 
489, L25-L28. 

Durech, J., Vokrouhlicky, D., Kaasalainen, M., Weissman, P., Lowry, S. C, 
Beshore, E., Higgins, D., Krugly, Y. N., Shevchenko, V. G., Gaftonyuk, 
N. M., Choi, Y.-J., Kowalski, R. A., Larson, S., Warner, B. D., Marshalk- 
ina, A. L., Ibrahimov, M. A., Molotov, I. E., Michalowski, T., Kitazato, K., 
Sep. 2008b. New photometric observations of asteroids (1862) Apollo and 



25 



(25143) Itokawa - an analysis of YORP effect. Astronomy & Astrophysics 
488, 345-350. 

Emery, J. P., Kelley, M. S., Fernandez, Y. R., Hergenrother, C. W., Crane, 
K. T., Ziffer, J., Campins, H., Lauretta, D. S., Drake, M. J., Oct. 2012. 
Thermal and Physical Characterization of the OSIRIS-REx Target Aster- 
oid (101955) 1999 RQ36. In: AAS/Division for Planetary Sciences Meet- 
ing Abstracts. Vol. 44 of AAS/Division for Planetary Sciences Meeting 
Abstracts, p. #102.05. 

Gehrels, T., Roemer, E., Taylor, R. C, Zellner, B. H., Mar. 1970. Minor 
planets and related objects. IV. Asteroid (1566) Icarus. The Astronomical 
Journal 75, 186-195. 

Giorgini, J. D., Benner, L. A. M., Ostro, S. J., Nolan, M. C, Busch, M. W., 
Jan. 2008. Predicting the Earth encounters of (99942) Apophis. Icarus 193, 
1-19. 

Giorgini, J. D., Ostro, S. J., Benner, L. A. M., Chodas, P. W., Chesley, S. R., 
Hudson, R. S., Nolan, M. C, Klemola, A. R., Standish, E. M., Jurgens, 
R. F., Rose, R., Chamberlin, A. B., Yeomans, D. K., Margot, J.-L., Apr. 
2002. Asteroid 1950 DA's Encounter with Earth in 2880: Physical Limits 
of Collision Probability Prediction. Science 296, 132-136. 

Hudson, R. S., Ostro, S. J., Jurgens, R. F., Rosema, K. D., Giorgini, J. D., 
Winkler, R., Rose, R., Choate, D., Cormier, R. A., Franck, C. R., Frye, 
R., Howard, D., Kelley, D., Littlefair, R., Slade, M. A., Benner, L. A. M., 
Thomas, M. L., Mitchell, D. L., Chodas, P. W., Yeomans, D. K., Scheeres, 
D. J., Palmer, P., Zaitsev, A., Koyama, Y., Nakamura, A., Harris, A. W., 
Meshkov, M. N., Nov. 2000. Radar Observations and Physical Model of 
Asteroid 6489 Golevka. Icarus 148, 37-51. 

Konopliv, A. S., Asmar, S. W., Folkner, W. M., Karatekin, 6., Nunes, D. C, 
Smrekar, S. E., Yoder, C. F., Zuber, M. T., Jan. 2011. Mars high resolu- 
tion gravity fields from MRO, Mars seasonal gravity, and other dynamical 
parameters. Icarus 211, 401-428. 

La Spina, A., Paolicchi, P., Kryszczyhska, A., Pravec, P., Mar. 2004. Retro- 
grade spins of near-Earth asteroids from the Yarkovsky effect. Nature 428, 
400-401. 



26 



Marsden, B. G., Sekanina, Z., Yeomans, D. K., Mar. 1973. Comets and 
nongravitational forces. V. The Astronomical Journal 78, 211. 

Micheli, M., Tholen, D. J., Elliott, G. T., May 2012. Detection of radiation 
pressure acting on 2009 BD. New Astronomy 17, 446-452. 

Mignard, F., 2003. Observations of Small Solar System Bodies with GAIA. 
In: IAU Joint Discussion. Vol. 19 of IAU Joint Discussion. 

Milani, A., Chesley, S. R., Sansaturio, M. E., Bernardi, F., Valsecchi, G. B., 
Arratia, O., Oct. 2009. Long term impact risk for (101955) 1999 RQ. Icarus 
203, 460-471. 

Moyer, T. D., 2003. Formulation for Observed and Computed Values of Deep 
Space Network Data Types for Navigation. Wiley-Interscience, Hoboken, 
NJ. 

Mueller, M., Aug. 2012. Surface Properties of Asteroids from Mid-Infrared 
Observations and Thermophysical Modeling. ArXiv e-prints. 

Nugent, C. R., Mainzer, A., Masiero, J., Grav, T., Bauer, J., Sep. 2012a. 
The Yarkovsky Drift's Influence on NEAs: Trends and Predictions with 
NEOWISE Measurements. The Astronomical Journal 144, 75. 

Nugent, C. R., Margot, J. L., Chesley, S. R., Vokrouhlicky, D., Aug. 2012b. 
Detection of Semimajor Axis Drifts in 54 Near-Earth Asteroids: New Mea- 
surements of the Yarkovsky Effect. The Astronomical Journal 144, 60. 

Pravec, P., Harris, A. W., Sep. 2007. Binary asteroid population. 1. Angular 
momentum content. Icarus 190, 250-259. 

Rozitis, B., Green, S. F., Jun. 2012. The influence of rough surface thermal- 
infrared beaming on the Yarkovsky and YORP effects. Monthly Notices of 
the Royal Astronomical Society 423, 367-388. 

Standish, E. M., 2000. Recommendation of DE405 for 2001 Mars Surveyor 
and for Cassini. Tech. Rep. IOM 312.F-00-107b, Jet Propulsion Laboratory. 

Standish, E. M., Hellings, R. W., Aug. 1989. A determination of the masses of 
Ceres, Pallas, and Vesta from their perturbations upon the orbit of Mars. 
Icarus 80, 326-333. 



27 



Taylor, P. A., Margot, J.-L., Vokrouhlicky, D., Scheeres, D. J., Pravec, 
P., Lowry, S. C, Fitzsimmons, A., Nolan, M. C, Ostro, S. J., Benner, 
L. A. M., Giorgini, J. D., Magri, C, Apr. 2007. Spin Rate of Asteroid 
(54509) 2000 PH5 Increasing Due to the YORP Effect. Science 316, 274-. 

Vokrouhlicky, D., Jul. 1998. Diurnal Yarkovsky effect as a source of mo- 
bility of meter-sized asteroidal fragments. I. Linear theory. Astronomy & 
Astrophysics 335, 1093-1100. 

Vokrouhlicky, D., Apr. 1999. A complete linear model for the Yarkovsky 
thermal force on spherical asteroid fragments. Astronomy & Astrophysics 
344, 362-366. 

Vokrouhlicky, D., Chesley, S. R., Matson, R. D., Jun. 2008. Orbital Identifi- 
cation for Asteroid 152563 (1992 Bf) Through the Yarkovsky Effect. The 
Astronomical Journal 135, 2336-2340. 

Vokrouhlicky, D., Milani, A., Chesley, S. R., Nov. 2000. Yarkovsky Effect on 
Small Near-Earth Asteroids: Mathematical Formulation and Examples. 
Icarus 148, 118-138. 

Will, C. M., Mar. 1993. Theory and Experiment in Gravitational Physics. 
Cambridge University Press. 

Yeomans, D. K., Antreasian, P. G., Barriot, J.-P., Chesley, S. R., Dunham, 
D. W., Farquhar, R. W., Giorgini, J. D., Helfrich, C. E., Konopliv, A. S., 
McAdams, J. V., Miller, J. K., Owen, W. M., Scheeres, D. J., Thomas, 
P. C, Veverka, J., Williams, B. G., Sep. 2000. Radio Science Results Dur- 
ing the NEAR-Shoemaker Spacecraft Rendezvous with Eros. Science 289, 
2085-2088. 



28 



Table 2: Scmimajor axis, eccentricity, absolute magnitude, physical and Yarkovsky pa r amete rs, and observational information 
for selected NEAs. The physical quantities for 1999 R Q36 are from Cheslev et al. ( 2012T). Goleyka's obliquity 7 is from 
Hudson et al.l (l2000j) . Apollo's from Durech et al. ( 2008bl) . Nyx's from Drummond and Wisniewski ( 1990l ). Ra-Shalom's from 
Durech et aiT (|2012T) . YORP's from iTavlor et al 



12007), Geographos' from iDurech et al.l (|2008al) . and Toro's from Durech 



(personal communication). Diameters D and rotation periods P are from the European Asteroid Research Node (EARN) 
database (http://earn.dlr.de/nea/). Bond albedo A was computed fro m geometric albedo py and slope parameter G (both 
from the EARN database when available) by A = (0.290 + 0.684G)pv (jBowell et al.l . Il989h . 
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HDP 
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A 2 


SNR S da/dt da/dt 


Diff. Observed Radar 




[au] 




[km] [hr] 






[lCT 15 au/d 2 ] 


[10~ 4 au/Myr] Nugent et al. 


a arc apparitions 



(101955) 1999 RQ 36 1.13 0.20 20.6 0.49 4.29 0.01 175° 

(152563) 1992 BF 0.91 0.27 19.6 0.51 NA 0.02 NA 

(6489) Golevka 2.52 0.60 19.1 0.27 6.03 0.23 137° 

2009 BD 1.01 0.04 28.2 NA NA NA NA 

(1862) Apollo 1.47 0.56 16.0 1.40 3.06 0.12 162° 

(2062) Aten 0.97 0.18 16.9 1.30 40.77 NA NA 

(3908) Nyx 1.93 0.46 17.3 1.00 4.43 NA 20° 

(2100) Ra-Shalom 0.83 0.44 16.1 2.24 19.80 0.05 162° 

(10302) 1989 ML 1.27 0.14 19.4 0.25 19 NA NA 

1999 MN 0.67 0.67 21.4 NA 5.49 NA NA 

(2340) Hathor 0.84 0.45 19.9 0.30 NA 0.24 NA 

(6037) 1988 EG 1.27 0.50 18.7 0.40 2.76 NA NA 

(37655) Illapa 1.48 0.75 17.8 NA 2.66 NA NA 

(85953) 1999 FK 2 i 0.74 0.70 18.0 0.59 NA NA NA 

(1685) Toro 1.37 0.44 14.3 3.00 10.20 0.05 144° 

2005 ES70 0.76 0.39 23.6 NA NA NA NA 

(54509) YORP 1.00 0.23 22.6 0.10 0.20 NA 178° 

(283457) 2001 MQ 3 2.23 0.46 18.9 NA NA NA NA 

(1620) Geographos 1.25 0.34 15.2 2.47 5.22 NA 153° 

(65679) 1989 UQ 0.91 0.26 19.3 0.73 7.73 NA NA 

(162004) 1991 VE 0.89 0.66 18.0 NA NA NA NA 



-45.49 ± 0.23 197.7 1.0 -18.99 ± 0.10 



-18.9 ± 0.2 0.4 1999-2012 



-24.01 ± 1.21 
-15.88 ± 1.52 



19.8 
10.4 



-1164.01 ± 138.76 8.3 



-3.34 ± 0.52 


6.1 


-15.41 ± 2.45 


6.3 


29.95 ± 5.39 


5.6 


-10.97 ± 2.25 


4.9 


90.48 ± 16.38 


4.6 


50.79 ± 11.39 


4.5 


-24.71 ± 5.66 


4.4 


-32.66 ± 8.19 


4.0 


-13.99 ± 3.65 


3.8 


-10.62 ± 2.33 


3.7 


-2.83 ± 0.77 


3.7 


-97.23 ± 29.27 


3.3 


-74.93 ± 23.68 


3.2 


-44.63 ± 14.33 
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-4.25 ± 1.40 


3.0 


-36.66 ± 12.23 
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18.30 ± 6.13 


3.0 



0.5 
0.2 

0.4 
0.2 
0.9 
1.3 
1.1 

1.0 
0.4 
0.3 
0.6 
0.6 
0.3 
0.1 

0.3 
0.3 
1.1 
0.5 
1.2 
0.7 



-11.55 ± 0.58 
-6.62 ± 0.64 

-493.39 ± 58.81 
-1.70 ± 0.26 
-6.29 ± 1.10 
11.61 ± 2.09 
-6.31 ± 1.30 

34.71 ± 6.28 
47.12 ± 10.56 
-14.33 ± 3.28 
-16.39 ± 4.11 
-11.27 ± 2.94 
-10.38 ± 2.28 

-1.27 ± 0.34 

-55.58 ± 16.73 
-33.60 ± 10.61 
-16.02 ± 5.14 
-1.82 ± 0.60 
-17.51 ± 5.84 
14.75 ± 4.94 



-12.84 ± 1 1.1 
-5.74 ± 0.7 0.9 



NA 
-2.3 ± 0.2 
-7.5 ± 2.4 
12.9 ± 2.7 
-5.4 ± 1.5 

35.3 ± 7.1 

NA 
-14.5 ± 3.5 
-14.2 ± 4.3 

NA 
-10.44 ± 1.5 
-1.4 ± 0.7 

NA 

-35.63 ± 10.5 
NA 
-2.5 ± 0.6 
NA 
NA 



NA 
1.8 
0.4 
0.4 
0.5 

0.0 

NA 
0.0 
0.4 
NA 
0.0 
0.2 

NA 
0.1 

NA 
0.8 
NA 
NA 



1953- 2011 
1991-2011 

2009-2011 
1930-2012 
1955-2012 
1980-2012 

1975- 2012 

1989-2012 

1999- 2010 

1976- 2012 
1988-2011 
1994-2012 
1971-2012 
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2005-2011 

2000- 2005 
1951-2011 
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1954- 2011 
1954-2012 



1999, 2005, 
2011 
NA 

1991, 1995, 
2003 
NA 

1980, 2005 
1995, 2012 
1988, 2004 
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2000, 2003 
NA 
NA 
NA 
1988 
2012 
NA 

1980, 1988 
2012 
NA 

2001, 2004 
NA 

1983, 1994 
NA 
NA 



CO 

o 



Table 3: Same as Table [2] fo r less reli a ble de tections. Icarus' obli quity is from Gehrels et al. (1970), Eros' from lYeomans et al.l 
( 2000h . Nereus' from Brozovic et al. ( 20091 ). and Betulia's from Drummond and Wisniewskil ( 199oT ). All the other physical 
quantities are from the EARN database. 
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Figure 9: Density as a function of thermal inertia for asteroids Golevka, Apollo, Ra- 
Shalom, Toro, YORP, and Geographos. The rectangles correspond t o reasonable value s 
of p according to the taxonomic type and to a reasonable range of r (|Delb6 et all . l2007h . 
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